Hadronic molecules for charmed and bottom baryons near 

thresholds 



Yasuhiro Yamaguchi 1 , Shunsuke Ohkoda 1 , Shigehiro Yasui 2 , and Atsushi Hosaka 1 
1 Research Center for Nuclear Physics (RCNP), 
Osaka University, Ibaraki, Osaka, 567-0047, Japan and 
2 KEK Theory Center, Institute of Particle and Nuclear Studies, 
High Energy Accelerator Research Organization, 
1-1, Oho, Ibaraki, 305-0801, Japan 
(Dated: January 31, 2013) 

Abstract 

We study hadronic molecules formed by a heavy meson and a nucleon, DN and D*N (BN 
and B*N) systems. Respecting the heavy quark symmetry and chiral symmetry, we consider the 
DN-D*N (BN-B* N) mixing induced by the one boson exchange potential including the tensor 
force. We find many bound and resonant states with J p = 1/2^, 3/2 ± , 5/2^ and 7/2 _ in isospin 
singlet channels, while only a few resonant states with J p = l/2~ in isospin triplet channels. The 
analysis of DN and D*N (BN and B*N) molecules will be useful to study mass spectra of excited 
charmed (bottom) baryons with large angular momenta, when their masses are close to the DN 
and D*N (BN and B*N) thresholds. 

PACS numbers: 12.39.Hg, 14.20.Lq, 14.20.Mr, 14.20.Pt 
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I. INTRODUCTION 



Recently, many exotic hadrons observed in experiments are attracting our interest. They 
are expected to have unexpected structures, because their properties, such as quantum num- 
bers, masses, and decay patterns, cannot be explained by the ordinary picture of hadrons, 
baryons as qqq and mesons as qq. As one of new aspects of these structures, the possibility 
of multiquarks or hadronic molecules have been extensively discussed. It is considered that, 
particularly near thresholds, hadronic molecules such as loosely bound or resonant states of 
meson-meson or meson-baryon emerge as one of hadronic configurations realized in nature. 
As typical candidates, the meson- meson molecules for A (3872) and have been studied in 



the charm and bottom quark regions 



The hadronic molecules as bound or resonant states formed by D or D* meson (B or B* 
meson) and a nucleon (N) are also interesting. Since the DN and D*N (BN and B*N) 
molecules contain ordinary flavor structures of three quarks, they are intimately related to 
the ordinary heavy baryons, A c , S c , S* and their excited states for the charm sector, and 



Ab, Sb, Sb an d their excited states for the 



studied in the quark model extensively 
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oottom sector M . These excited states have been 



12j. However, near the DN and D*N (BN and 



B*N) thresholds, it is expected that the wave functions of the quark model states have a 
large overlap with the component of DN and D*N (BN and B*N), and the properties of 
such baryon states are strongly affected by DN and D*N (BN and B*N) states. There 
would be even such states that are dominated by the molecular components. 

In the DN and D*N (BN and B*N) systems, there are two important symmetries; the 
heavy quark symmetry and chiral symmetry. The heavy quark symmetry manifests in mass 
degeneracy of heavy pseudoscalar mesons P = D, B and vector mesons P* = D*, B* in the 
heavy quark mass limit flfl. Indeed, the mass splitting between a D (B) meson and a 
D* (B*) meson is small; 140 MeV for D and D* mesons and 46 MeV for B and B* mesons. 
The small mass splittings in the heavy flavor sectors should be compared with the large 
mass splittings, ~ 400 MeV for K and K* mesons, and ~ 600 MeV for tt and p mesons, in 
the light flavor sectors. Because of this, the heavy quark symmetry introduces the mixing 
of PN and P*N states, where both P and P* mesons are considered on the same footing 



as fundamental degrees of freedom in the dynamics 15|. Then, physical states are given as 



a superposition of PN and P*N states. For convenience, in the following, we denote P^ 



to stand for P or P*. Thanks to the mixing of states in p(*>N systems, a long range force 
between a P^ meson and a nucleon N is supplied by one pion exchange with the same 
coupling strengths for PP*ir and P*P*tt vertices 16H18] . 

It is interesting that the pion exchange can exist in p(*)jV systems. As a matter of fact, 
it is known that the pion exchange is important for the binding of atomic nuclei 19j. There 
the tensor force of the pion exchange that mixes channels with different angular momenta, 
i.e. L and L ± 2, yields a strong attraction to generate a rich structure of bound as well as 
resonant states. Because the pion exchange is a long range force, it becomes more significant 
for states with larger orbital angular momenta. In fact, it was shown that, in the nucleon- 
nucleon scattering, the phase shifts in the channels with large orbital angular momenta are 



reproduced almost only by the pion exchange potential [20|, |2l| . The importance of the pion 
exchange was also discussed in NN systems. It was investigated that the properties of bound 
and resonant states in NN systems are dramatically changed, if the tensor force is switched 



off 
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23J. In our previous studies, it was also shown that the pion exchange potential 



lfHial. DD 



played a significant role for hadronic molecules such as DN and BN baryons 
(BB) mesons (e.g. Zb meson) 24| and DD (BB) mesons 25|. 

In reality, however, the hadronic molecules of P^N do not necessarily correspond to the 
observed states, because there should be couplings not only to three quark states but also 
to other meson-baryon states such as 7rA c , 7rE c and 7rE* (vrAb, vrSb and vrE^), and so on. 
However, we may expect that such couplings are small for DN and BN baryons near the 
thresholds. The reasons are that the wave functions of hadronic molecules are spatially large 
as compared to the conventional three quark states, and that the transitions, e.g. from DN 
(BN) to 7rS c (7rSb), are suppressed by a heavy quark exchange. From those points of view, 
in the present discussion, we focus on the role of the D^N and B^N sectors and study the 
bound or resonant states generated by D^N and B^N. Recently, there have been several 
studies of baryon states as D&N (B^N) and £>WE C (B^E h ) 
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molecules. 



The present study is different from them in that we focus on the states near the thresholds 
and with a large angular momentum. 

This paper is organized as follows. In Sec. HI1 we briefly describe the interaction between 

a 

£)(*) or #(*) 

meson and a nucleon. In Sec. Illfl we solve the coupled channel Schrodinger 
equations numerically for bound and resonant states. We discuss the results of the D^N 
and B^N molecule in Sec. IIVt and summarize the present work in the final section. 



II. INTERACTIONS 



Following our previous studies 116 
tween and N. The interaction Lagrangians 
heavy quark symmetry and chiral symmetry [39 



18) , we consider n, p and u exchange potentials be- 
br heavy meson vertices are given by the 



4Q|, 



Cvhh = ~i/3Tr [H b v^) ba H a ] + iATr [H b a^F^(p) ba H a ] 



(1) 
(2) 



where the subscripts tt and t> are for the pion and vector (p and u) meson interactions. The 
heavy meson field H is defined by H a = [Pa^l^ ~ -Pa 75] with the four-velocity of a 
heavy meson, where the subscript a is for light flavors, u, d. The conjugate field is given by 
H a = 70-^70. The axial current of light flavors is written by = | — 0^^) with 

£ = exp (iTr/fir) and the pion decay constant f n = 132 MeV. The pion field is defined by 



7T 



/ 


7T° 










7T~ 


V 





(3) 



and the vector meson field by 



.9v . 
V2 



P 



(4) 
(5) 



P "71/ 



where gy ~ 5.8 corresponds to the gauge coupling constant from hidden local symmetry (41]. 
The vector meson field tensor is written by F^ u (p) = d^pu — d u p^ + [p^, p v ]. The coupling 
constants g n , (3 and A are summarized in Table [I] These coupling constants are essentially 



the same as those given in Refs. 



17 



18], except for the signs of g n , (3 and A for vertices of 
tt and to, which are reversed due to G-parity transformation between Z)W (£?(*)) and 



(B^**). In Eqs. ([Q) and (EJ, we consider the static approximation 



:i,o). 



For vertices of w, p and 00 mesons to a nucleon, we employ the Bonn model (42] as 

V2ig nNN N l5 frN , (6) 



C 



vNN 



V2g 



vNN 



K 



N^N+- Na^Vp^N 

2m 



N 



(7) 
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TABLE I. Masses and coupling constants of mesons a = ir, p, co in Refs. [171 . [18]. 



Meson 


m a 


9-K 


P 


A [GeV" 1 ] 


91nn/^ 


K 


TT 


137.27 


-0.59 






13.6 




P 


769.9 




0.9 


0.59 


0.84 


6.1 




781.94 




-0.9 


-0.59 


20.0 


0.0 



where N = (p, n) T is the nucleon field. The coupling constants for nucleons are summarized 
in Table HI 

When we consider the ir, p and u exchange potentials between meson and nucleon, 
we take into account internal structure of hadrons by introducing the monopole-type form 
factors at each vertex; 



F a (A,q\ 



A 2 - m 2 
A 2 + |<f 2 | 



(8) 



with the mass m a , three momentum q of the incoming meson a (= ir, p, u), and the cutoff 
parameter A. Here, we introduce two cutoff parameters for a heavy meson (£>M and #W) 
and a nucleon. The cutoff parameter A^r for the nucleon is determined to reproduce the 
properties (the binding energy, the scattering length and effective range) of the deuteron. 
The cutoff parameter Ap for heavy meson = D^*\ is determined by the ratios of 
matter radii of the heavy meson and nucleon where the ratio is obtained by a quark model 

rical values are summarized in Table HTl 



calculation 



In the present study, we employ two types of potential: the tt exchange and the ixpuj 
exchanges. By comparing the results from them, we show the importance of the one pion 
exchange potential. 



TABLE II. Cutoff parameters of a nuc 
for i?M meson) as employed in Ref. [17 



ucleon (Ajv) and heavy mesons (A^> for DM meson and Ap 



Potential A N [MeV] A D [MeV] Ap [MeV] 



7T 

irpuj 



830 
846 



1121 
1142 



1070 
1091 
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TABLE III. Various coupled channels for a given quantum number J p . 



J p 


channels 


1/2- 


PN( 2 S 1/2 ) 


P*N( 2 S 1/2 ) 


P*N( 4 D l/2 ) 




1/2+ 


PN( 2 P 1/2 ) 


P*N{ 2 P l/2 ) 


P*N( 4 P 1/2 ) 




3/2" 


PN( 2 D 3/2 ) 


P*A( 4 S 3/2 ) 


P*N( 4 D 3/2 ) 


P*N( 2 D 3/2 ) 


3/2+ 


PN{ 2 P V2 ) 


P*N( 2 P 3/2 ) 


P*N( 4 P 3/2 ) 


P*N( 4 F 3/2 ) 


5/2" 


PN( 2 D n/2 ) 


P*N( 2 D 5/2 ) 


P*N( 4 D 5/2 ) 


P*N( 4 G 5/2 ) 


5/2+ 


PN( 2 F 5/2 ) 


P*N( 4 P 5/2 ) 


P*N( 2 F 5/2 ) 


P*N( 4 F 5/2 ) 


7/2- 


PN( 2 G 7/2 ) 


P*N( 4 D 7/2 ) 


P*N( 2 G 7/2 ) 


P*N( 4 G 7/2 ) 


7/2+ 


PN( 2 F 7/2 ) 


P*N( 2 F 7/2 ) 


P*N( 4 F 7/2 ) 


P*N( 4 H 7/2 ) 



III. NUMERICAL RESULTS 

We consider the states with J p = 1/2+, 3/2+, 5/2+ and 7/2+ (total angular momentum 
J and parity P) with isospin 1 = and 1. Various states with J p are expanded by coupled 
channels of 2S+1 Lj (spin S and orbital angular momentum L) as summarized in Table. IHIl 
The meson exchange potentials among various channels and kinetic terms are summarized 
in Appendix |A] Using the hamiltonian composed of the kinetic and potential terms, we solve 
the coupled-channel Schrodinger equations for PN and P*N channels numerically. 

A. Isospin singlet (/ = 0) 

We present the results for the isosinglet state (1 = 0). The energies of bound and resonant 
states are summarized in Table HVl and also presented in Figs. [1] and [2j The results are shown 
for the 7r and irpu potentials, and the corresponding energy levels are connected by arrows 
in Figs. [1] and [2j The numerical values are measured from the DN and BN thresholds, 
respectively. 

First we discuss the states with J p = 1/2+ and 3/2", and then those with J p = 3/2+, 
5/2+ and 7/2+. In fact, it will turn out that the results in those two categories have different 
behaviors. 
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1. J p = l/2 ± and 3/2" 



For (J, J p ) = (0, 1/2"), we find bound states in both of D^N and For D^N, the 

binding energies are —14.4 MeV for the n potential and —82.5 MeV for the npu potential. 
The relative radii are 1.51 fm and 0.86 fm, respectively. As expected, for a larger binding 
energy the system becomes smaller. The tensor force from the ir exchange which is the main 
driving force of the DN-D*N mixing plays an important role to form the bound states. In 
fact, we have verified that neither bound nor resonant state exists when the tensor force 
from the ir exchange is switched off and the DN-D*N mixing is small. For J p = 1/2", 
the results in the ir and irpu potentials are very different as indicated in Figs. [1] and [21 
Since both p and u exchanges are attractive for D^N (B^N) system, the vector meson 
exchanges contribute to form the deeply bound state of the binding energy around 80 MeV. 
This contrasts with the previous result for the D^N and B^N systems of truly exotic 



channels, where t 
two potentials 117 



r e p and u exchanges play a minor role due to the cancellation of these 



We estimate the mixing ratio of various channels in the bound states as summarized in 
Table IVl We observe that, for J p = 1/2" states, the most dominant component is DN( 2 Si/ 2 ) 
with a fraction 71.8 %. The second dominant component is D*N{ 4l Di/ 2 ) with a fraction 20.4 
%. Therefore, the tensor force which mixes the 5*-wave in DN( 2 Si/2) and the D-w&ve in 
D*N( A Di/ 2 ) is important to provide a strong attraction. 

For B^N of (/, J p ) = (0, l/2~), we also obtain bound states. The binding energies are 
—57.8 MeV for the ir potential and —145.9 MeV for the irpuj potential, and the relative 
radii are 0.99 fm and 0.76 fm, respectively. Again the results of the two potentials are very 
different with the same reason for the D^N system. As compared to the D^N system, 
the binding energy in the B^N system is much larger, because heavier particles suppress 
kinetic energy. For the mixing ratios, we also find a similar tendency as we discussed for the 
D^N system; 56.1 % for BN( 2 S 1/2 ), 13.3 % for B*N( 2 S l/2 ) and 30.6 % for B*N( 4 D 1/2 ). 

Let us move to (/, J p ) = (0, l/2 + ) state. For D^N, we find one resonance near the DN 
threshold when the 7r potential is used. The resonance energy is 1.4 MeV and the half decay 
width is 0.2 MeV as shown in Table IIV1 In the scattering state, we define the resonance 



energy E Te by an inflection point of the phase shift 43(, and the decay width T is given by 



T = 2 / \d5 1 'dE) E=E rc - When the npu potential is used, we find a bound state with a binding 



TABLE IV. Properties of bound and resonant states of D^N and B^N systems. The energies 
E are either pure real for bound states or complex for resonant states. The complex energies for 
resonances are written as E re — iT/2 where E re is the resonance energy and T/2 is the half-width. 
The binding and resonance energies are measured from the lowest threshold (DN and BN). Root 
mean square radii are shown only for bound states. 



I(J P ) 


Potential 


DN 


BN 


E [MeV] (r 2 ) 1 / 2 [f m ] 


E [MeV] (r 2 ) 1 / 2 [f m ] 


0(1/2") 


TT 
TTpOJ 


-14.4 1.51 
-82.5 0.86 


-57.8 0.99 
-145.9 0.76 


0(1/2+) 


TT 

irpu) 


1.4 -iO. 2 
-81.5 0.85 


-83.8 0.92 
-185.0 0.75 


0(3/2") 


TT 


63.5 - £7.9 

-13.7 0.89 


-38.7 0.99 
-127.8 0.76 


0(3/2+) 


TT 

irpu) 


23.8 - £118.1 
26.0 - i44.2 


12.9 - £15.5 

-2.6 1.81 


0(5/2-) 


TT 
TTpOJ 


153.6 - £671.9 
160.0 - £375.4 


63.7 -£177.6 
71.3 - £102.8 


0(5/2+) 


TT 

irpu 


160.8 - £3.1 
137.0 - £7.6 


46.0 -£1.1 
20.0 -£0.2 


0(7/2") 


TT 
TTpOJ 


217.7 -£182.4 
220.8 -£109.1 


85.6 - £74.5 
87.5 - £46.7 


0(7/2+) 


TT 
TTpOJ 


no 
no 


no 
no 


1(1/2-) 


TT 
TTpOJ 


no 

147.2 - £105.5 


no 

50.7- £75.5 



energy —81.5 MeV and with a relative radius 0.85 fm. The mixing ratios of the bound state 
are 38.8 % for DN( 2 P 1/2 ), 6.0 % for D*N( 2 P 1/2 ) and 55.2 % for D*N( 4 P 1/2 ), as shown in 
Table. [V] Interestingly, D*N( 4 Pi/ 2 ) is the most dominant channel, although the mass of 



E 



220.8 - ilOQ.l 



7/2" 



160.0 - i375.4 



D*N 
(2949 MeV) 



J?/?: 



217.7 -£182.4 
7/2" 

153.6 - 4671.9 
5/2- 



160.8-13.1 



137.0-i7.6 



'5/2+ 



DN 
(2807 MeV) 



63.5 - «7.9 



"3/2" 



-13.7 




-14.4 



1/2" 



26.0 - i44.2 23.8 -£118.1 
372^T4~gr_3/2+ 



^Ti - 1/2 + 



7TpW 7T 
P = - 



TVp UJ TV 
P = + 



FIG. 1. Energies of bound and resonant states of D^N for 1 = with positive parity (P = +) 
and negative parity (P = —). The energies are measured from the DN threshold. The results for ir 
and irpuj potentials are shown, where corresponding states are connected by arrows. The binding 
energies are given as real negative value, and the resonance energies E re and decay widths T are 
given as E re — iT/2. 

D*N is heavier than the mass of DN. This is because the attraction of tensor force is the 
strongest for D*N{ A P\/2) channel. 

For B&N of (/, J p ) = (0, 1 /2+), we find bound states for both cases when the tc and npu 
potentials are used. The binding energies are —83.8 MeV for the tt potential and —185.0 
MeV for the irpu potential, and the relative radii are 0.92 fm and 0.75 fm, respectively. 
For the npu potential, the mixing ratios of the bound state are 28.4 % for BN( 2 Pi/ 2 ), 7.7 
% for B*N( 2 P 1/2 ) and 63.9 % for B*N( i P 1/2 ). In this case again, B*N( 4 P 1/2 ) is the most 
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E 



B*N 
(6265 MeV) 

BN 
(6219 MeV) 



87.5-i46.7 7/2 85 . 6 _ i74 . 5 



71.3-il02.8 



5/2" 



63.7 - (177.6 



-145.9 



.7/2" 
'5/2- 




-38.7 3 / 2 ~ 
^S - 1/2- 



46.0 -il.l 



20.0 - £0.2 5 /' 2+ 
i</2+ 12.9 - ?15.5 

3/2+ 

~ ^2~6 3/2+ 



-185.0 




-83.; 



1/2+ 



p = - 



■KpUJ 7T 
P = + 



FIG. 2. Energies of bound and resonant states of B^N for 1 = with positive parity (P = +) and 
negative parity (P = — ). The energies are measured from the 5iV threshold. The same convention 
is used as Fig[TJ 

dominant component regardless of its heavy mass. 

For (J, J p ) = (0, 3/2 _ ), we find a resonant DN state for the 7r potential with the res- 
onance energy 63.5 MeV and the half decay width 7.9 MeV as shown in Table IIVI In 
contrast, for the irpu potential, we find a bound DN state with the binding energy —13.7 
MeV and the relative radius 0.89 fm. For the bound state, the mixing ratios are 19.8 % 
for DN( 2 D 3/2 ), 62.8 % for D*N( 4 S 3/2 ), 14.2 % for D*N( 4 D 3/2 ) and 3.2 % for D*N( 2 D 3/2 ). 
Thus, D*N( A S 3 / 2 ) is the dominant channel although its mass is heavier than the mass of 
DN( 2 D 3 / 2 ); once again a large attraction due to the tensor force is provided. 

For B^N state of (I,J P ) = (0,3/2~), we obtain bound states for both cases the ir 
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TABLE V. Mixing ratio of each channel in the bound D^N and B^N states for J p = 1/2+ and 
3/2+ with 1 = when the irpu potential is employed. 



1/2- 


PN( 2 S 1/2 ) P*N( 2 S 1/2 ) P*N( i D 1/2 ) 


1/2+ 


PN( 2 P 1/2 ) P*N( 2 P 1/2 ) P*iV( 4 P 1/2 ) 


DN 
BN 


71.8% 7.8% 20.4% 
56.1% 13.3% 30.6% 


DN 
BN 


38.8% 6.0% 55.2% 
28.4% 7.7% 63.9% 


3/2 - 


PN( 2 D 3/2 ) P*N( 4 S 3/2 ) P*N( 4 D 3/2 ) P*N( 2 D 3/2 ) 


3/2+ 


PN( 2 P 3/2 ) P*N{ 2 P 3/2 ) P*7V( 4 P 3/2 ) P*N( A F 3/2 ) 


DN 
BN 


19.8% 62.8% 14.2% 3.2% 
14.6% 64.7% 16.7% 4.0% 


DN 
BN 


71.2% 6.7% 7.5% 14.6% 



and irpuj potentials are used. The binding energies are —38.7 MeV for the it potential and 
— 127.8 MeV for the irpuo potential with relative radii 0.99 fm and 0.76 fm, respectively. 

2. J p = 3/2+, 5/2+ and 7/2+ 

For (I,J P ) = (0,3/2+), we find resonant states above the threshold. For D^N, the 
resonance energies are 23.8 MeV for the tt potential and 26.0 MeV for the irpuj potential. 
The half decay widths are 118.1 MeV and 44.2 MeV, respectively. As compared to the cases 
of 1/2+ and 3/2~, the results of the 7r and irpu potentials are not very much different, as 
indicated in Fig. [TJ Since the wave functions are extended due to larger orbital angular 
momenta of P-wave (L = 1) and F-wave (L = 3) for J p = 3/2+, the long range potential 
of the 7T exchange dominates, while the short range potentials from p and u exchanges 
are suppressed. For B~(*'N, when the tt potential is used, we find a resonant state whose 
resonance energy is 12.9 MeV and half decay width is 15.5 MeV, while when the irpui 
potential is used, we find a loosely bound state of a binding energy —2.6 MeV and a relative 
radius 1.81 fm. The mixing ratios are 71.2 % for BN( 2 P 3/2 ), 6.7 % for B*N( 2 P 3/2 ), 7.5 % 
for B*N( 4 P 3/2 ) and 14.6 % for B*N( 4 F 3/2 ). The most dominant component is BN( 2 P 3/2 ), 
and the second dominant one is B*N( 4 F 3 / 2 ). 

For (/, J p ) = (0, 5/2"), we obtain resonant states for both cases when the 7r and irpu 
potentials are used. For D^N, the resonance energies are 153.6 MeV for the tt potential and 
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160.0 MeV for the Tipu potential, which are above the D*N threshold. The corresponding 
half decay widths are 671.9 MeV and 375.4 MeV, respectively. The difference between the 
results of the ir and npoo potentials is once again small, due to the same reason as before 
with large angular momenta. For B^N state, we also find resonant states above the B*N 
threshold. The resonance positions are 63.7 MeV for the n potential and 71.3 MeV for 
the npuj potential. The corresponding half decay widths are 177.6 MeV and 102.8 MeV, 
respectively. 

For (J, J p ) = (0, 5/2 + ), we find resonant states with narrow widths for both cases when 
the Ti and ixpoj potentials are used. For D^N, when the 7r potential is used, the resonance 
energy is 160.8 MeV and the half decay width T/2 = 3.1 MeV. When the irpu potential 
is used, the resonance energy is 137.0 MeV and the half decay width V/2 — 7.6 MeV. For 
B&N, we also find resonances whose energies are 46.0 MeV for the it potential and 20.0 
MeV for the 7rpu potential, with the corresponding half decay widths 1.1 MeV and 0.2 MeV, 
respectively. Again, the results for the tt and npu potentials are similar. 

For (J, J p ) = (0, 7/2"), we obtain resonances above the D*N and B*N thresholds. For 
D&N, there exist resonances at 217.7 MeV for the 7r potential and at 220.8 MeV for the 
Tvpoo potential, with half decay widths 182.4 MeV and 109.1 MeV, respectively. For B^N, 
there also exist resonances whose energies are 85.6 MeV for the it potential and 87.5 MeV 
for the npu potential, with half decay widths 74.5 MeV and 46.7 MeV, respectively. 

Finally, for (I, J p ) = (0, 7/2 + ), we find no bound nor resonant state. 
B. Isospin triplet (I = 1) 

Let us move to the states of isospin triplet (1 = 1). We summarize the results for D^N 
and for B^N in Table HVl and show the energy levels in Fig. [3j As a result, we find resonant 
states only for J p = 1/2" when the Trptu potential is employed. For D^N, the resonance 
energy is 147.2 MeV and the half decay width is 105.5 MeV. For B^*'N state, we obtain the 
resonance whose energy is 50.7 MeV and half width 75.5 MeV. The reason that there are 
not many states in isospin triplet channel can be understood as follows; as compared to the 
isosinglet channel, the attractive force in the isotriplet channel is weak due to small isospin 
factor; r P ■ r N = —3 for isosinglet and r P ■ r N = 1 for isotriplet. 
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E 



147.2 - il05.5 



D*N 
(2949 MeV) 



Ll/2- 



DN 
(2807 MeV) 



B*N 
(6265 MeV) 

BN 
(6219 MeV) 



50.7 - £75.5 



1/2: 



P 



P = + 



P 



FIG. 3. Energies of bound and resonant states of D^N and 2?WjV for 7 = 1 when the npoj 
potential is used. The same convention is used as Fig{TJ 

IV. DISCUSSION 

A. Contribution from short range interactions 

By solving coupled channel Schrodinger equations for PN and P*N channels, we find 
many bound and resonant states for 1 = and few resonances for 1=1. For these states, 
the tensor force of the n exchange potential plays a significant role to produce them. When 
we ignore the P*N channels and solve the Schrodinger equation only with PN channels, we 
find no bound state nor resonance. Thus, the PN-P*N mixing and the rr exchange potential 
play an important role to generate a rich structure of molecular states. 

The importance of the ir exchange potential stands out more in the bottom sector. The 
small mass difference between B and B* mesons helps to yield the strong attraction because 
it induces the strong BN-B*N mixing with the tensor force. Furthermore, the B^*'N has 
heavier reduced mass and hence the kinetic term is suppressed. For these reasons, the 
binding energies of B&N states are larger than these of D^N states. 

When we compare the results of the rr potential and of the npu potential, they are quite 
different for J p = l/2 ± , 3/2~ with 1 = 0, where the p and u exchange potentials become 
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important to produce attraction. For D^*>N and B^N states, both p and u exchange 
potentials are attractive, and hence they increase the binding energy significantly. On the 
other hand, for J p = 3/2 + , 5/2 ± and 7/2", the results for the irpu potential are similar 
to those for the ir potential. For large J states, the it exchange potential plays a dominant 
role to generate bound and resonant states, while the p and u exchanges play only a minor 
one. This is attributed to the fact that these states have large orbital angular momenta. If 
relevant channels include large orbital angular momenta, the wave functions tend to extend 
spatially, and the long range force, namely the 7r exchange potential, becomes important 
while the short range force is suppressed. 



B. Emergent hadronic molecules 

In the present study, we discuss the molecular structure formed by the P^N bound and 
resonant states. However, the hadronic molecular picture is not applicable to deeply bound 
states with small radii. For such states, the constituent hadrons, namely (B^) meson 
and nucleon, overlap each other, and therefore we have to consider short range effects such 
as an internal structure of hadrons, channel couplings to conventional three quark states 
and so on. As a naive criterion for the hadronic molecule, we have shown the relative radii 
of the bound states as discussed in Ref 



251 ]. If the size of the bound state is larger than 
twice of typical hadron size (namely 1 fm), the state could be well described by a molecular 
structure. For resonant states, we identify the states as the hadronic molecule. According to 
the criterion, for the up to potential, only the bound state for J p = 3/2 + of B^N constructs 
a hadronic molecule, where the relative radius is 1.81 fm. Contrary to this, the bound states 
for J p = 1/2^ and 3/2~ with 1 = which have a small radius and a large binding energy 
are not described as simple molecules. We need to consider the short range effects including 
various channel couplings to do more realistic discussions. 



V. SUMMARY 

We have investigated heavy baryons as hadronic molecules formed by a heavy meson and 
a nucleon, D^N and B^N. The interaction is given by the meson exchange potential with 
respecting the heavy quark and chiral symmetries. Many bound states and resonances are 
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found in isosinglet channel both for D^*'N and B^N systems. In contrast, there are few 
resonances in isotriplet channel due to small isospin factor. For J p = l/2 ± and 3/2~, we 
have found deeply bound states far below DN or BN threshold with a small radius. In these 
states, the vector meson exchange potential yields a strong attraction. In order to perform 
more realistic discussions for such compact states, we need to consider further effects from 
internal structures of constituent hadrons and channel couplings to quark intrinsic states. 
On the other hand, for J p = 3/2 + , 5/2 ± and 7/2~ with large orbital angular momenta, 
the one pion exchange potential is favored and we have obtained loosely bound states and 
resonances near the thresholds. There the short range interactions are rather inactive and 
the long range force of the one pion exchange dominates, where the tensor force plays an 
important role to generate these states, and they compose the hadronic molecular structure. 

It is expected that the hadronic molecular states from a heavy meson and a nucleon are 



studied in many accelerator faci 



annihilation process in FAIR 



ities; the D meson production in the antiproton-nuc 



44l | and the heavy ion collision in RHIC and LHC 



45 



con 



46). 



Furthermore, the exotic baryons are investigated in J-PARC in the coming future. 
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Appendix A: Potentials and kinetic terms 



Interaction potentials are derived by using the Lagrangians in Eqs. ([I]), ([2]), §6§ and ([7]) 



as shown in Ref. 
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. The potentials for the coupled channel systems are given in the 



matrix form of 3 x 3 for J p = 1/2^ and of 4 x 4 for the other states. 
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V3V£ -2Vg 

1 



\/lVn 



_J_T/ 7r _ 



T 



where 



v c 



9n9nNN 1^ _ 



The fp and are the isospin matrices for and iV. The functions C, 
T m = T(r; m) are given by 



C(r; m) 
T(r; m)S 12 {f) 



d 3 q 



m 



(27r) 3 q 2 + m 2 



e^F(A P ,q)F(A N ,q), 



(A8) 



d 3 g 



^i 2 (g)e^F(Ap,g)F(A^,g; 



(A9) 
C(r; m) and 

(A10) 
(All) 



(27r) 3 <f 2 + m 2 

with Si2(:c) = 3((?i • x)(<?2 • x) — <?i • <?2, and F(A, q) denotes the form factor given in Eq. (jBJ). 
The vector meson exchange potentials (u = p, u) are 



/ 



V l/2~ 



V 1 " 



2y/3Vg V6V£ 
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Tj 



(A13) 
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/ 



l/ 7/2^ 



v c 



2y/3Vg VS'-iVS 
1 



V3 T 



5 v? 



V% Vg' + 2V%--V2 



c 



V£' + 2V£ + -Vv 



(A19) 



where Vq', Vq and V£ are defined as 



Vg = — jL „ C mo T P ■ T N 



V p 
v c 



9v9 p nnH 1 + K ) 1 ri - - 



V2m 



N 



r /P _ 9v9 p nnK 1 + «0 !^ 



v c 



y/2m 

9v9ioNnP 

V2ml 



Jp ■ tn 



N 



Vf7 — P= TT^n 



V, 



T 



y/2m,N u 
9v9ujNN^ 1 
\f2~Tn j\f 3 



(A20) 

(A21) 
(A22) 
(A23) 
(A24) 
(A25) 



Finally, the kinetic terms are given by 

1 



if 1/2- = diag 



A , — 

2m p 2m p* 



1 A 1 

A +Ampp*, 



2m p* 



A 2 +Ampp« 



(A26) 



K 1/2 + = diag 



2m, 



■Ai, 



2m P * 



Ai +Ampp*. 



2m P * 



Ai +Ampp. 



(A27) 



A' 



3/2- 



diag 



2m, 



■A 2 , 



1 



2m p* 



Aq +Ampp». 



1 



2m p* 



A 2 +Ampp* 



2m p» 



A 2 +Ampp* 



(A2£ 



^3/2+ = diag f -— — Ai, -77^— Ai +Amp P ,,-— Ai +Am PP , , 
\ 2m p 2m p* 2m p* 

A 3 +Ampp* J , (A29) 



2mp 
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K 5/2 - = diag ( -__A 2 , -77^— A 2 +Am P p.,--3- A 2 +Am PP , , 
\ Imp Imp* Imp* 

A 4 +Am P p. J , (A30) 



2m p 



K 5/2 + = diag ( — — ^A 3 , -— — A x +Am PP ., -— - — A 3 +Am PP . , 
\ Imp Imp* Imp* 

A 3 +Am Ppi ) , (A31) 



2m p 



K 7/2 - = diag ( -— 5— A 4 , -— ?— A 2 + Am PP ,,--^- A 4 +Am PP , , 
\ 2m p 2m p* 2m p* 

A 4 +Am PP . ) , (A32) 



2fhp 

K 7/2 + = diag ( -— _A 3 , -77-^— A 3 +Am PP .,--A A 3 +Am PP . , 
\ 2m p 2mp* 2mp* 

A 5 +Am PP *) , (A33) 



2fhp 

where A; = d 2 / dr 2 + {2/ r)d / dr — /(/ + l)/r 2 , m P (,j = m^mpi*) /(m^+mpi*)), and Am PP « 
mp* — mp. The total Hamiltonian is then given by Hjp = Kjp + Vjp. 
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